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Abstract 

We consider the optimal trade execution strategies for a large portfolio of single 



stocks proposed by Almgren ( 2003 1 . This framework accounts for a nonlinear impact 



of trades on average market prices. The results of Almgren (20031 are based on the 



assumption that no shares of assets per unit of time are trade at the beginning of 
the period. We propose a general solution method that accomodates the case of a 
positive stock of assets in the initial period. Our findings are twofold. First of all, 
we show that the problem admits a solution with no trading in the opening period 
only if additional parametric restrictions are imposed. Second, with positive asset 
holdings in the initial period, the optimal execution time depends on trading activity 
at the beginning of the planning period. 
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1 Introduction 



The execution of large trades in financial markets requires the balance between risks and 
costs. The main risk concerns the lack of availability of a counterparty, which can lead to 
a delay in the execution of a transaction. In order to guarantee a fast trade execution, a 



trader may incur additional costs. As clarified by Hasbruock and Schwartz (19881, a trader 
faces a choice between a 'passive' and an 'actice' execution strategy. 

Given this background, the available models of optimal execution assume that the 
trading activity of individual investors has an impact on the average price prevailing in the 
market. The transaction costs are characterized by parametric forms that replicate stylized 
facts documented in the market microstructure literature (e.g. see Kraus and Stoll, 1972). 



Almgren and Chriss (1999, 2000) and Konishi and Makimoto (2001) provide examples 
of optimal strategies for the execution problem in the stock market. Their models assume 
that the transaction cost per share is a linear function of the number of shares of assets 
traded. The only source of uncertianty consists in the volatility of the stock price. 



Almgren (2003) suggests that the linearity assumption is largely at odds with reality. 



First, the average liquidity premium on stocks tends to be either a convex or a concave 
function of the traded size. This depends on the counterparty's perception about the 
reason for the trade, namely on whether it is driven by liquidity or information needs (see 



Huang and Stoll, 1997). Moreover, the liquidity premium is related to the risk of finding 



a counterparty. In other words, the lower the probability of finding a counterparty in the 
market, the higher the liquidity premium. 



In this paper, we review the optimal transaction strategy proposed by Almgren (2003). 



We show that the solution method used by Almgren (2003) is ill-posed. The reason is 
that it is based on the assumption that no shares per unit of time are exchanged at the 
beginning of the period. We use an approach based on the Gauss hypergeometric function 
to solve for the case of positive initial trades. Our results differ strongly from those of 



Almgren (2003). First of all, the problem admits a solution with no trading in the opening 



period only if additional parametric restrictions are imposed. Second, with positive initial 
trading, the optimal execution time depends on trading activity in the initial period. 

This note is organized as follows. Section 2 provides a selected discussion of the 
literature on optimal trade execution. Section 3 outlines the structure of the problem. 
Section 4 proposes a general solution method for positive initial values of the velocity. 
Section 5 concludes. Finally, in Appendix A, we discuss the general method for the solution 
of second order differential equations with a Gauss hypergeometric function. 
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2 A selected overview of the literature 



Rebalancing portfolios of assets requires executing trades in the marketplace. With the 
advent of algorithmic trading and access availability to many alternative trading venues, 
investors have dedicated increasing resources to the scheduling of trades. The practical 
setup of the problem is rather intuitive. An investor has a target number of, say shares 
that it intends to sell or buy within a given time frame. The decision problem requires to 
compute how many shares to place or demand in the market at each point in time within 
the trade horizon. The aim of the investor is to minimize the execution costs. These are 
typically measured as the difference between the price obtained from the market and a 
benchmark price for the transaction. 

There are multiple relevant dimensions to the execution problem. Several contributions 
have showed that the liquidity premium is time- varying. The reason is that it is determined 
by the availability of traders willing to act as counterparties, namely traders willing to buy 
or sell a given quantity of an asset at a desired price. However, as the presence of traders 
willing to 'take the other side' of a trade is uncertain, any trading is characterized by 
execution risk. 

Another relevant aspect is related to the fact that market illiquidity generates 
transaction costs. This typically takes the form of a large spread between bid and ask 



prices (Huang and Stoll 1997). Therefore, as noticed by Wagner and Banks (1992), the 
minization of transaction costs is a key aspect of the portfolio optimization problem. As 



documented in various studies including Chakravarty (2001), Holthausen, Leftwich and 



Mayers (1990) and Kraus and Stoll (1972), large trades do impact market prices and, thus 



affecting the bid-ask spreads. 

Asset price volatility is a source for execution risk. The reason is that it affects the 
probability of finding a suitable counterparty. Hence, it affects the sucessfulness of a trading 
strategy execution. The recent literature has focused on the specific aspect of volatility, 
namely the increased uncertainty in execution price incurred by rapid execution of large 
share blocks. In fact, Hasbruock and Seppi (2001) show that liquidity fluctuates due to 
intrinsic variations in market activity independently of trade size. 

Based on the considerations outlined earlier, what are the properties of an 'optimal' 



execution strategy? What defines an 'optimal execution price' ? Bertsimas and Lo ( 1998 1 
argue that 'best' execution can be thought of as a dynamic strategy that minimizes 
'liquidation' costs. They show that dynamic programming techniques can be used fruitfully 
to characterize these strategies. 

Almgren and Chriss (2000) compute optimal trajectories for trading prices that are 
obtained by balancing market impact costs. The optimal profiles provide a motivation for 
low execution speedj^ This results arises from the balance between the need to reduce 

^The literature has proposed two main alternative bechmarks. These consist of average prices that 
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the expected value of execution, and the need to minimize the adverse effects of market 
volatihty. While the first factor provides a reason for slower execution, the second factor 
lays the ground for rapid execution. That would, in fact, reduce execution risk in the 
form of the variance of execution cost with respect to the benchmark price. In short, early 
execution reduces execution risk, whereas a delayed execution is more geared towards 
minimizing execution costs. Evidently the degree of investor risk aversion determines how 
early within the trading horizon the execution starts. The shape of the schedule instead 
depends on the form of the assumed market impact model. 



Konishi and Makimoto (2001 1 makes the assumption that the market impact of trading 
is a linear function of trade size. In this case, the optimality frontier representing 
the combinations of minimized costs and market volatility has an analytical solution. 
Value-at-risk utility funcations are then used to select the first-best solution. This choice 
of utility function provides a natural testing ground for the concept of liquidity-adjusted 
VaR, which explicity considers the best trade-off between volatility risk and liquidation 
cost. Almgren (2003) generalizes the results of Konishi and Makimoto (2001) to the case 
of nonlinear functions for the market impact of trades. In this framework, the assumption 
is that the market impact cost per share follows a power law function of the trading rate. 

3 The optimal execution problem 

We follow the general framework of Almgren and Chriss (2000). At time i = 0, an investor 
holds X shares of an asset. The problem is to sell these shares by by time t = T. We 
should stress that this is the statement of a general framework. In fact, the initial size X 
can either be positive or negative. In the first case, the investor needs to schedule a selling 
program. In the second case, the investor looks at a buying program. In this paper, for 
simplicity, we focus on the case X > 0. 

The execution problem consists in minizing the market impact of trades subject to both 
initial and terminal conditions. In mathematical terms, the model proposed by [Almgren" 



(2003) delivers the following optimization problem: 



( r 

min / x) dt 



(3.1) 

^ x(0) = X, x{T) = 
where F{x, y) is the market impact function of trades: 

F(x,y) = -7xy + 7?(-y)'=+i + AaV, 7, r/. A, /o 0. (3.2) 



materialize within the trading horizon, and are characterized as a time-weighted and volume-weighted 
average price. 
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The problem is to determine the optimal function x{e) so as to minimize a chosen cost 



functional. Using Beltrami identity, (see Kamien and Scwartz 1991 section 5, page 31): 



F{x,x) — x——{x,x) = constant 
dy 



(3.3) 



evaluating the constant of integration at the end time T we are lead to the differential 
equation 



Acj^ — kri{—x) 



k+l 



-kr]{—x{T)) 



k+l 



(3.4) 



Almgren (2003) proposes a solution based only on the "elementary" case vq = 0. In this 
note, we show how to compute the constant Vq"*"^ {—x{T))''~^^ using the initial condition 
x(0) = X. 



3.1 The case k = l 

The case /c = 1 is straightforward since it gives rise to a linear ordinary differential equation, 
whose solution is better found starting from the classical Euler Lagrange equation: 



d dL dL 



> 2Xa^x -2r]x = 
imposing the boundary conditions x{0) = X, x(T) = we find: 



(3.5) 



x{t) = X 



sinh(^a(r-t) 
sinh(^ar 



(3.6) 



It is worth noting that for k = 1 (3.6 ) gives a minimizer of (3.1 ) since in this case Legendre 
condition reads: 



(x, x) = 2ry > 



(3.7) 



Turning back to equation (3.4) for general k, we write it solving with respect to x : 



fc+l A(j2 2\'=+i 

[x(0)=X 
In this general case observe that, since: 



(3.8) 



dy- 



-{x, x) = r]k{k + 1)(— x) 



k-l 



(3.9) 
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since solution to (3.8) is decreasing we infer, for the Legendre condition, the minimaUty of 
the estremal x{t). 



4 A general solution method for the case of positive initial 
trades 



In this section, we propose a solution method that holds when there are positive stock 



trades in the initial period, namely for vq > 0. This solution x to (3.8) is implicitely 
defined by 



X 



,fe+i 



ki] 



1 

' fc+i 



dz = t. 



(4.1) 



Integral in the right hand side of (4.1) can be evaluated by means of the Gauss 



hypergeometric function 2F1 whose definition and basic properties are given in the 



appendix. After some changes of variables which allows to rewrite (4.1) as: 



1 
2^ 



X 



_ 1 
s 2 



1 , Xa2x2 2 



ds — X 

k+T Jo 



_ 1 

S 2 



1 I Ao-2a;2 2 



— ds 



(4.2) 



fe+i 



we can finally use the integral representation Theorem, see equation (A. 4) in the appendix, 
for the hypergeometric function to obtain 



Vo 



X2F1 



1 _i_ 

2' k+1 

3 
2 



X2F1 



1 _j_ 

2' k+1 

3 
2 



Xa'^x 



2^2 



k+1 



(4.3) 



To obtain x{t) from equation (4.3) observe that the function 

Xa^x^ \ 



X I—)- X 2F1 



1 _j_ 

2' fc+1 



(4.4) 



is strictly decreasing function for values of the independent variable > and so is possible 



to revert and obtain x{t) from (4.3), if needed, numerically or, better, using the Lagrange 
power series reversion when: 



Xa^x^ 



kr]VQ~^^ 



< 1 



(4.5) 
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But in (4.3) there is no determination of vq which is essential for the full solution of the 
problem. If we limit to assign some specific values for vo as in the "easy case" t;o = we 
lose control on the initial value x{0) = X. 

The way to obtain the final velocity in order to fit with the initial value x(0) = X is 
explained below. We use the so called "shooting method" as presented for instance in |Stoer 



and Burlish (1993) section 7.3.1 pages 502-507. The starting point is the Euler Lagrange 



equation with initial values in T 



r]k{k + 1 
x{T) = 0, x{T) 



x[—x 



(4.6) 



To solve (4.6) we use the change of variables u := x, y := x following Murphy (1960) 
section 2.3 pages 160-161 and, since 

Ay 



x = y 



(4.7) 



equation (4.6) is transformed in 
2Aa2 



' dy ^ 
dn 

2/(0) 



r]k{k + 1) 

-Vq 



i-y) u 



Since (4.8) is separable we can integrate it, so that, returning back to the original variables 
we find: 



,A:+1 



H -X' 

ijk 



1 

fc+i 



(4.9) 



[x{T) = 



Equation (4.9) is separable and the relevant integration needs again the hypergeometric 



integral. Discarding some tedious computations we find out that solution to (4.6) is defined 
implicitely by: 



X 



Vo 



2-r 1 




Xa'^x 



2^2 



kr]VQ 



k+l 



T 



(4.10) 



In order to meet the second initial condition x(0) = X we see that vo must satisfy: 



X 



2^ 1 




Xa^X^ 



kvvo 



k+l 



T 



(4.11) 
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Note that, being assigned all parameters A, o", X, rj, k equation (4.11 ) is an equation in the 



sole unknown vq. Of course such equation has to be treated numerically: once the value 



of Vq is detected it has te be inserted in (4.3) to obtain solution to (3.1). 

As an (easy) example take X = a = X = ri = T= l and k = 1/2. In the plot below 



we represent the left hand side of (4.11 ) and the value of T at the right hand side. Using 
Mathematica® we find numerically vq = 0.671525. 

T 



1.5 
1.0 
0.5 



Figure 1: Plot of equation (4.11) 



Afterwards we put this value in (4.3) and we plot the relevant function x{t) 




0.2 0.4 0.6 0.8 1.0 



Figure 2: Plot of solution to (3.1) 
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4.1 The Almgren zero-speed case 



When Almgren takes vq = the initial value problem (4.6) has to be considered with zero 
initial conditions x{T) = x{T) = 0. Observe that the quadrature formula arising from 



(4.9), which is equivalent to (4.6) reads in this case as: 
dz 







t 



+ jjfc 



dz 



T 







(4.12) 



fc+l 

rjk 



Zk+l 



but this means that the convergence condition < 1 has to be imposed. Moreover since 
we assume k to be positive this means that the zero speed assumption is well posed only if 



k > 1: for k < 1 there are not solutions of the optimization problem (3.1 ) with zero speed 



Moreover if we evaluate the integral at the left hand side of (4.12) for k > 1 we find: 



k + l 
k-1 



kr] 



1 

k-1 
Xk+1 



T-t 



x{t) 



(k 



k + l 

l){T-t)\ 



k + l 



krj J 



(4.13) 
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so that at t = we have: 
x(0) = 



k + l 



(4.14) 



this means that we are not free to consider an arbitrary value of x{0) having assigned the 
speed in i = T. 

4.2 The case k = l 



In the case k = 1 we have provided two solutions of (3.1): the first follows from the 



straightforward integration of the linear case, see equation (3.6), while the second stems 



from the hypergeometric implicit solution expressed by equation (4.3). Of course the two 



solution to (3.1) are, as matter of fact, the same. This can be understood recalling the 



following property of 2F1, (see Abramowit and Stegun, 1964, entry 15.1.7, page 556): 

(4.15) 



2-r 1 







/ 2' 2 




I i 


-) 







In (z + Vl + z'^ 



arcsinh z 



In such a way, in this particular case we can use this identity to solve (4.10) with respect 
to t obtaining: 



xit) 



v,sinh({T-t),/^ 



(4.16) 



To compare this hypergeometric solution with (3.6) we have to evaluate vq from (4.11): 
^/XaX 



Vo 



(4.17) 



and substitute in (4.16) getting: 



x{t) = X 



(4.18) 



which is noting else but (3.6) 
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5 Conclusion 



Rebalancing large portfolios of stocks requires taking into account two peculiar issues. The 
first one is related to the market impact of trades, which generates transaction costs. The 
second issue arises from the risk of finding counterparties willing to trade at the desired 
price. Both empirical and theoretical considerations suggest that the market impact of 



trades is typically nonlinear. Almgren (2003) proposes an optimal execution strategy 
that minimizes the tradeoff between volatility risk and transaction costs while taking into 
account this form of nonlinearity. 



In this paper, we review the optimal liquidation strategy of Almgren (2003). We show 



that the solution method used by Almgren (2003) is ill-posed. The reason is that it is 
based on the assumption that no shares per unit of time are traded at the beginning of the 
period. We use an approach based on the Gaussian hypergeometric function to solve for 



the case of positive initial trades. Our results differ strongly from those of Almgren (2003). 
First of all, the problem admits a solution with no trading in the opening period only if 
additional parametric restrictions are imposed. Second, with positive initial trading, the 
optimal execution time depends on trading activity in the initial period. 
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A Solution of second order differential equations with a 
Gauss hypergeometric function 

The linear second order differential equation in the unknown u = u{t) 



t{l-t)u+[c - {a + b + l)t] u-abu = 



(A.l) 



is known as Gauss hypergeometric equation. Parameters a, b, c are not functions of the 
independent variable t and can be in general complex number. Searching for a power series 



solution of (A.l) it can be seen that 

" (a)„(6)„ 



2F1 



a, b 



n=0 



ni 



(A.2) 



where we use the Pochhamer symbol {x)n, € N defined as: 
{x)o := 1 

{x)n = x(x + l)(x + 2) • • • (x + n - 1) 



(A.3) 



is the solution of (A.l) such that n(0) = 1, tt(0) = ab/c. Power series defining 2F1 
converges for |t| < 1 and to continue the hypegeometric function 2F1 it is useful the integral 
representation ascribed to Leonhard Euler but really due to Adrien Marie Legendr^ 



a, b 
c 



r(c) 



r(c-a)r(a) 



1 s"-i(l - s)=-"-i 



{1-ts)'' 



-ds, 



(A.4) 



where Rec > Re a > 0, \t\ < 1, and the Euler-Legendre integral (Gamma function) is 
defined for x > by: 



r(x) 



(A.5) 



A proof of the integral representation theorem and a good presentation of the Gauss 
hypergeometric function can be found at Seaborn ( 1991[ ), the integral representation 
theorem is treated at section 10.7, pages 184-185, formula (10.39). Integral representation 
theorem provides an extension to the region where the complex hypergeometric function 
is defined, namely for its analytical continuation, to the (almost) whole complex plane 
excluding the half-straight line (l,oo). This function was first introduced in dynamical 



economics in a generalization of the Solow Swan model due to Mingari Scarpello and Ritelli 



(2003), while Boucekkine and Ruiz-Tamarit (2008) use it in the Lucas-Uzawa model. 



^Exercices de calcul integral, II, quatrieme part, section 2, Paris, 1811. 
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